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The present study describes the genotypic distribution of 
rotaviruses (RVs) in an Indian bovine population with 
unexpectedly higher proportions of G3 alone or in 
combination of G8/G10. PCR-genotyping confirmed that 
39.4% (13/33) of the prevalent RVs were the G3 type while 
60.6% (20/33) were dual G3G10 or G3G8 types. P typing 
revealed that 93.9% (31/33) of the samples were P[ll] while 
6.1% (2/33) possessed a dual P[1]P[11] type. Sequence 
analysis of the VP7 gene from G3 strains viz. B-46, 0970, and 
BR-133 showed that these strains had sequence identities of 
90.5% to 100% with other bovine G3 strains. The highest 
identity (98.9% to 100%) was observed with RUBV3 bovine 
G3 strains from eastern India. The G3 strains (B-46, 0970, 
and BR-133) showed 97.5% to 98.8% sequence homologies 
with the Indian equine RV strain Erv-80. Phylogenetic 
analysis demonstrated that G3 strains clustered with bovine 
RUBV3 and J-63, and equine Erv-80 G3. Overall, these 
results confirmed that the incidence of infection by RVs with 
the G3 genotype and mixed genotypes in the bovine 
population was higher than previously predicted. This 
finding reinforces the importance of constantly monitoring 
circulating viral strains with the G3 genotype in future 
surveillance studies. 

Keywords: bovine group A rotavirus, distribution, G and P type, 
India, PCR genotyping 



Introduction 

Rotavirus (RV) is one of the leading causes of diarrhea in 
neonatal animals, humans, and birds. Infection with RVs 
has resulted in colossal economic losses in the livestock 
industry [8]. These losses are much greater in developing 
countries like India where socio-economic conditions 
favor the survival of virus in the environment [11]. RVs are 
icosahedral non-enveloped viruses belonging to the family 
Reoviridae with a genome comprising of 1 1 segments of 
double-stranded (ds) RNA, each of which encodes one of 
the six structural (VP1 ~4, 6, and 7) or six non-structural 
(NSP1-6) proteins [7]. The complete infectious RV virion 
consists of three concentric layers: the inner core, inner 
capsid, and outer capsid [7]. Three out of the six structural 
proteins (VP4, VP6, and VP7) have major antigenic 
determinants and are used for RV classification into groups 
(A~G), sub-groups (SGI/II), and genotypes/serotypes 
(G&P) [7]. Group A rotavirus (RVA) is the most 
predominant serogroup present in bovines and other 
mammalian species throughout the world [7]. 

Based on neutralizing antibodies against the outer capsid 
proteins VP7 and VP4, two independent glycoprotein (G) 
and protease sensitivity (P) serotype systems have been 
established. However, sequence-based genotyping 
[4, 1 1 , 1 2, 1 4, 1 5,3 0] targeting the VP7 (G type) and VP4 (P 
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type) protein genes has become more a common method 
for RVA typing due to difficulties in generating serotyping 
reagents and performing the time-consuming assays. To 
date, 25 G (Gl ~25) types and 32 P (P[l ~ 32]) types have 
been identified for different RV strains [3,6,19,32]. 
Genotypes in RVs circulating among different animal 
populations have been examined in order to design 
effective vaccines and assess RV cross-species infectivity. 
Changing patterns in prevalence of bovine G and P types 
have been observed over time. The most predominant 
genotypes in bovines are G6, G8, and G10 with sporadic 
incidences of Gl, G3, G5, or G15. The most prevalent P 
types are P[l], P[5], and P[l 1] with sporadic occurrence of 
P[14], P[15], P[17], and P[21] [2,7,9,11,17,21,23,25,26]. 

In India, G10 (50~60%) has been the predominant 
bovine RV genotype over the last 1 0 years followed by G6 
(25-35%) [11,14,20,22,27,28,34]. Varshney et al. [33] 
were the first to report on the existence of G3 RV in buffalo 
calves ( 1 0.2%) in central India. This was the second most 
predominant strain after G10. In 2007, Ghosh et al. [11] 
also reported a 3% incidence of G3 RVs in claves from 
eastern India, indicating that either the number of bovines 
infected with RV G3 genotypes are increasing or these 
animals have not been taken into account while conducting 
genotyping surveys in the country [14,20,22,27,28,34]. 
The present study describes G and P genotypes prevalent in 
Indian bovine populations from diverse tropical, sub- 
tropical/semi-arid, and temperate western Himalayan 
regions during 2007 ~ 20 10. Sequence and phylogenetic 
analyses of the G3 genotype presently circulating in the 
country were also performed. 

Materials and Methods 

Sample collection and preparation 

A total of 378 fecal samples were collected from cattle 
and buffalo calves up to 3 months old that were suffering 
from diarrhea. These animals belonged to organized and 
unorganized dairy farms in different ecoregions of India 
including temperate western Himalaya (Tarai and foothills 
of Himalaya, Uttarakhand state), tropical, semi-tropical, 
and semi-arid regions (Haryana, Uttar Pradesh, and 
Madhya Pradesh) of India. A 1 0% fecal suspension (w/v) 
was prepared with phosphate buffered saline (0.01 M, pH 
7.4; Sigma, USA). Samples were centrifuged at 2,000 x g 
for 10 min to remove coarse particulate matter and the 
upper aqueous layer was filtered through a 0.22-um pore 
filter (MDI, India) into a fresh tube. The suspension was 
stored at — 20°C until further use. Reference NCDV 
(G6P[1]) and UK (G6P[5]) cell culture-adapted bovine 
RVA strains were obtained from the National Research 
Centre on Equine (India). RNA extracted from these 
samples was used to optimize reverse-transcription 
(RT)-PCR and genotyping conditions. 



Nucleic acid extraction from the fecal suspensions 

Total RNA was extracted from 500 uL of fecal suspension 
using an equal volume of TriReagent-LS (Sigma, USA) 
following the manufacturer's instructions. RNA was 
dissolved in nuclease-free water (NFW) (Life 
Technologies, USA) in a final volume of 25 uL. RNA 
quantity and quality was assessed using a Nanodrop 
Spectrophotometer (ND-1000; Thermo Scientific, USA) 
and stored at — 20°C until further use. 

Viral RNA electrophoresis (RNA-PAGE) 

The extracted viral RNA was analyzed by PAGE 
performed as previously described by Laemmli [18] by 
using 5% stacking and 10% resolving gels. Briefly, 
extracted RNA (200 — 500 ng) was mixed with an equal 
amount of 2x RNA loading buffer (Thermo Fisher 
Scientific, USA) and loaded into the wells of the gels. 
Electrophoresis was performed at 120 V in 1 x Tris-glycine 
buffer (25 mM Tris, 192 mM glycine, pH 8.3). The gel was 
subsequently stained with silver nitrate as previously 
described by Svensson et al. [29]. All the samples were 
initially screened with PAGE. 

RT-PCR 

RT of the viral RNA was performed using a random 
hexamer primer (0.2 ug/uL; Thermo Fisher Scientific, 
USA). Initially, 50 ~ 100 ng of viral RNA, 0.5 uL random 
hexamer primer, and 2 \\L dimethyl sulphoxide (DMSO) 
were combined in a thin-walled 0.5-mL PCR tube 
(GenAxy, India). The reaction mixture was incubated at 
95°C for 5 min to denature the dsRNA strands and 
immediately snap chilled on ice. Next, 4 ixL of 5x RT 
buffer (Promega, USA), 2 ixL of 10 mM dNTPs (Thermo 
Fisher Scientific, USA), 40 U of RNase inhibitor (Life 
Technologies, USA), and 200 U of Moloney murine 
leukemia virus reverse transcriptase (Promega, USA) were 
added and the mixture was incubated at 37°C in a 
thermocycler for 1 h (HR-PCR-96G; Haier, China). The 
reverse transcriptase was denatured at 80°C for 5 min at the 
end of the incubation step. 

The first-round of PCR amplification of the VP7 gene 
was performed with 3 uL cDNA, 5 uL 10x PCR buffer, 2 
mM MgCl 2 , 2 uL DMSO, 1 uL dNTPs, 50 pmol of forward 
and reverse primers (Table 1), and 1.25 U of Taq DNA 
polymerase (Thermo Fisher Scientific, USA). The total 
reaction volume was brought to 50 uL with NFW. The 
first-round of amplification of full-length VP7 (1062 bp) 
consisted of an initial denaturation at 94°C for 4 min 
followed by 35 cycles of 1 min at 94°C, 2 min at 46°C, and 
2 min at 72°C, and final incubation at 72°C for 1 0 min. A 
similar PCR protocol was performed to amplify the partial- 
length VP4 gene (864 bp) for the first-round with the 
exception that the annealing temperature was set at 48°C. 
The PCR products were separated by conventional agarose 
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Table 1. Primers used for amplifying VP7 full-length and VP4 partial-length amplicons, and G and P typing 



Primer Sequence Position . „ _ Reference 

size (bp) 



VP7 full-length 


VP7-F 


5'GGCTTTAAAAGAGAGAATTTCCGTCTGG 3' 


1 


~28 


1062 


Taniguchi et al. [31] 




VP7-R 


5'GGTCACATCATACAATTCTAATCTAAG 3' 


1062 


-1036 






G typing 


G3 


5 CGTTTGAAGAAGTTGCAACAG 3' 


689 


-709 


373 


Gouvea a/. [12] 




G5 


5 GACGTAACAACGAGTACATG 3' 


760 


-779 


302 


Gouvea et al. [13] 




G6 


5 GATTCTACACAGGAACTAG 3' 


481 


-499 


581 


Gouvea et al. [13] 




G8 


5'GTCACACCATTTGTAAATTC 3' 


178 


-197 


884 


Gouvea et al. [12] 




G10 


5 ' ATGTCAGACTACARATACTGG 3' 


666 


-686 


396 


Iturriza-Gomara et al. [16] 


VP4 partial-length Bov4Com5 


5'TTCATTATTGGGACGATTCACA 3' 


1067 


-1088 


864- 






Bov4Com3 


5 CAACCGCAGCTGATATATCATC 3' 


1930 


-1909 




Isegawaeia/. [15] 


P typing 


P[l] 


5'TTAAATTCATCTCTTAGTTCTC 3' 


1526 


-1505 


460 






P[5] 


5 GGCCGCATCGGATAAAGAGTCC 3' 


1725 


-1704 


659 — 






P[ll] 


5'TGCCTCATAATATTGTTGGTCT 3' 


1398 


-1377 


332 





G (VP7) typing 



P (VP4) typing 



NA 
5' , 



Amplification products 
1 062 bp | 



First round amplicon 



884 bp |- 



Genotype G8 



581 bp|- 



Genotype G6 



396 bp |- 



Genotype G10 



*> 



/ / ^ 

NT i f f f £ 



Amplification products 

First round amplicon 
864 bp | 



o9> 



2359 
-I* 



332 bp |- 
460 bp |- 



H P[H] 
— I PHI 



^ Full-length primers 
■ Typing primers 



659 bp |- 



H P[5] 



Genotype G3 
373 bp | — 1 

Genotype G5 
302 bp | 1 

Fig. 1. Diagrammatic illustration of VP7 (G)- and VP4 (P)-specific primer positions and type-specific amplicons. 



gel electrophoresis (1%, w/v) at 100 V for 1 h in lx TAE 
buffer with 0.5 |ig/mL ethidium bromide, viewed with a 
UV transilluminator, and documented using Syngene gel 
documentation system (Syngene, UK). 

Genomapping 

Multiplex nested PCR genotype determination was 
performed as per the method and primers used by Gouvea 
etal. [12] forG3 (aET3)andG8 (aAT8), Gouvea et al. [13] 
for G5 and G6, and Iturriza-Gomara et al. [16] for G10. 
Primers used for the first-round of amplification and 
subsequent nested PCR for genotyping both the VP4 and 
VP7 genes are presented in Table 1. Genotyping was 



performed using 3 uL of the amplicons (1 : 100 dilution) 
from the first-round of PCR amplification and the typing 
primers (G3, 5, 6, 8, and 10) with the reverse primer used 
for generating the first-round VP7 gene product. 
Amplicons from the first-round of PCR (1 : 100 dilution) 
and the typing primers (P[l], [5], and [11]) with the 
forward primer (Bov4Com5) used for generating the 
first-round VP4 gene product were used for P typing. A 
diagrammatic sketch with specific amplicon size for the G 
and P types is presented in Fig. 1 . PCR cycling conditions 
for the second-round were same as those of the first-round. 
The amplicons were analyzed as described above. 
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Sequencing and sequence analysis 

Selected bovine RVA isolates (n = 3) with the G3 genotype 
from cattle in different geographical areas were amplified as 
described above. The amplified PCR products were excised 
from the gel and purified with a High Pure PCR Product 
Purification kit (Roche Diagnostics, Germany) with final 
elution in 20 ixL of elution buffer. The eluted PCR products 
were quantified in Nanodrop Spectrophotometer and 
sequenced with an automated sequencer (ABI 3130, USA) 
at the Foot and Mouth Disease Typing Facility Lab, Project 
Directorate on Foot and Mouth Disease, Indian Veterinary 
Research Institute, Mukteswar (India). Sequence data for 
these samples were entered into the NCBI GenBank 
database under the following accession numbers: B-46 
(HM235510), BR-133 (JF 720877), and 0970 (HQ 199897). 

The sequence chromatogram was visualized with 
sequencing analysis software (ABI, USA). MegaBLAST 
(NCBI, USA) was used to analyze the predicted sequence 
within the non-redundant nucleotide database and confirm 
the presence of the gene specific for RVs. RotaC web-based 
RV classification software (Katholieke Universiteit Leuven, 
Belgium) was also used for sequence confirmation. RVA 
sequences from animals in different geographical locations 
within India and the rest of the world were downloaded 
from the NCBI nucleotide database. Sequence analysis 
was performed using the LaserGene software package 
(DNASTAR, USA). Percent identity and divergence of the 
nucleotide and predicted amino acid sequences were 



determined to calculate the percentage homology with 
other G3 RV sequences from different geographical 
locations. Aphylogenetic tree was constructed using Mega 
4.0 software [30]. 

Results 

Electropherotyping 

PAGE analyses detected the RV genome in 33 out 378 
diarrheic fecal samples from bovine calves that were up to 
3 months old. All PAGE-positive RV samples showed a 
typical 4:2:3:2 "long" electropherogram migration 
pattern, characteristic of RVA. In the majority of the 
samples, genome segments 7, 8, and 9 co-migrated as a 
triplet. PAGE-positive samples were also confirmed by 
using RT-PCR targeting the group A-specific VP6 gene. 

RT-PCR and genotyping 

First-round amplification with generic primers of VP7 
and VP4 genes produced full-length 1062-bp amplicons 
for the VP7 gene and partial-length 864-bp products for the 
VP4 gene in all 33 samples. G typing produced amplicons 
373, 581, 884, and 396 bps in size, indicating the presence 
of G3, G6, G8, and G10 genotypes, respectively (Fig. 2A). 
None of the samples were found to have the G5 or G6 type 
RV. P typing of the samples produced amplicons 334, 453, 
and 661 bp in size, corresponding to P[l 1], P[l], and P[5] 
types, respectively (Fig. 2B). However, P[5] was not 



A 

M1 1 



2 3 4 5 6 7 8 9 10 11 12 13 M2 



1000 bp 



500 bp 



Fig. 2. Agarose gel electrophoresis showing amplicon sizes for the different G and P types. (A) G types, Lane Ml: 1-kb DNA ladder, 
Lane M2: 1 00-bp DNA ladder, Lane 1 : VP7 full-length amplicon ( 1 062 bp), Lane 2: genotype G6 (NCDV strain), Lanes 3 ~ 6 and 1 1 ~ 
13: G3 types, Lanes 7 and 8: G3G10 types (single sample showing amplification of two viz. G3 and G10 genotypes), Lanes 9 and 10: 
G3G8 types. (B) P types, Lane M: 100-bp DNA ladder, Lane 1 : P[5] (UK strain), Lane 2: P[l] (NCDV strain), Lanes 3, 4, 6 and 7: P[l 1] 
type, Lanes 5 and 8: P[1]P[11] dual P types. 
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Table 2. Distribution of G and P types in different ecoregions of India 



State 
(Ecoregion) 


Rotavirus 
positive 
samples 




O types (to j 


r types (to) 




G3 G6 


G8 Glfl Mixed 

V_J O V_J 1 \J 1VJ.1A.CU 


prn pr5i pnn 

r L x J r L J J r L 11 J 


Mi y pH 

IVllAvU 


Uttarakhand 


17 


9 (52.9) 


6G3G8 (35.3) 


16 


1P[1]P[11] 


(Temperate western Himalaya) 






2G3G10(11.8) 






Uttar Pradesh 


3 


1 (33.3) 


1 G3G8 (33.3) 


3 




(Tropical subtropical) 






1G3G10(33.3) 






Haryana 


4 




2 G3G8 (50) 


4 




(Tropical subtropical/semi-arid) 






2 G3G10 (50) 






Madhya Pradesh 


9 


3 (33.3) 


6G3G10(66.7) 


8 


1P[1]P[11] 


(Tropical subtropical/semi-arid) 












Total 


33 


13 (39.4) 


20 (60.6) 


31 (93.9) 


2(6.1) 








9/20 G3G8 (45) 












11/20 G3G10(55) 







detected in any of the field isolates (Fig. 2B). A number of 
samples were confirmed to possess dual G types, either 
G3G8 or G3G10. The distribution of different combination 
of genotypes found throughout the study is summarized in 
Table 2. Typing primers (G&P) also confirmed the 
reference bovine RVA as G6P[ 1 ] for the NCDV strain and 
G6P[5] for the UK strain (Figs. 2A and B). G3 was found 
to be the main G type circulating among the ecoregions of 
India we examined between 2007 -2010 with both single 
and/or dual infections (G3G8 or G3G10) observed (Fig. 
2A). Overall, 13 of the 33 typed samples (39.4%) were 
identified as G3, and 20 (60.6%) were identified as mixed 
(G3G8 or G3G10) G types (Table 2). Additionally, 31 
(93.9%) of the 33 P typed samples had the P[l 1] genotype 
and two (6.1%) exhibited dual P[1]P[11] types (Table 2). 
The majority of the mixed G types were G3G10 (55%, 
11/20) following G3G8 (45%, 9/20) types. 

Distribution of G and P types 

The distribution of the G types showed that 52.9% were 
G3 and 47% were mixed G types. Among these, 35.3% 
were G3G8 (6/17) and 11.8% G3G10 (2/17) in the 
temperate region of western Himalaya while 25% of G3 
and 75% mixed G types were found in the tropical, 
subtropical, and semi-arid regions. Among the mixed 
types, the most commonly detected G combination was 
G3G10 (55%). P types found in the temperate region of 
western Himalaya were mostly P[ 1 1 ] (94. 1 %) while 5.9% 
were the dual P type (P[l] P[H])- In the tropical, 
sub-tropical, and semi-arid regions (Table 2), 93.8% 
(15/16) were P[ll] type and 11.1% (1/16) were the dual P 
type (P[1]P[11]). 

Sequence analysis 

The multiplex nested PCR typing results were further 



confirmed with sequencing, BLAST, and automated 
genotyping RotaC software analysis of selected isolates. 
As shown in Fig. 3, VP7 coding region sequence analysis 
of the bovine G3 isolates (B-46, 0970 and BR- 133) showed 
that the predicted amino acid sequence identities with other 
bovine G3 RVA strains were between 90.5% (Bo/UK/1984 
strain UKg9P GQ225779) and 98.8% (Bo/Ind/2006/ 
RUBV3 EF200549). Identity was lower with bovine 
strains from the UK (90.5% to 92%) while two other 
bovine G3 strains isolates previously obtained from two 
different parts of India showed higher identities of 96.3% 
(J63 strain from central India) to 100% (RUBV3 strain 
from eastern India). Sequence identities of these three G3 
isolates were between 89.9% and 95.1% with feline and 
canine G3 isolates, 93.3% and 94.5% with simian isolates, 
and 91.7% and 93.1% with equine isolates from Greece, 
Japan, and Germany. Equine G3 isolates from India had 
higher identities of 93.9% to 98.8%. All the G3 bovine 
isolates showed the lowest identity with porcine strains 
(82.8% to 91.4%). Sequence identities of these bovine 
strains with human G3 isolates from different parts of the 
world were between 89% and 94.5% (Fig. 3). 
UK bovine G3 RVs showed a maximum of 1 83 variations 
at the nucleotide level corresponding to changes in the 
amino acid sequence at 26 positions compared to B-46, 
0970, and BR-133 G3 strains. The RUBV3 strain from 
eastern India was more conserved as only one change was 
noticed at nucleotide level while there was no change at 
amino acid level. The J-63 G3 strain from central India 
showed 20 changes at the nucleotide level and eight 
changes at amino acid level (positions 22, 75, 146, 147, 
211, 213, 221, and 242). Compared to the J-63 strain, the 
B-46 G3 strain isolated from the same location showed 
only 1 0 variations at the nucleotide level and four changes 
(59, 71, 76, and 3 17) at amino acid level. 
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B-46 



Percent identity 
0970 



0.15 



0.10 



0.05 



0.00 



BR133 





NA 


AA 


NA 


AA 


NA 


AA 


Hu l nailana zuu4 uMn^ AYAU/vyz 


84.1 


93.6 


86.2 


94.8 


85.8 


94.5 


uan Italy zuuo Kvoz/yb Arz/iuyu 


84.3 


93.9 


86.2 


95.1 


85.9 


94.8 


tq inaia zuu / trv i uo uijyts 1 4/ y 


84.9 


93.9 


86.1 


95.1 


85.7 


94.8 


tq japan zuub riu/o AbU4b4b4 


83 


91.4 


84.7 


92.6 


84.5 


92.3 


bq UGrmany zuuo 4/bbUo AY/ouyz4 


81.8 


91.7 


83.3 


92.9 


83.2 


92.6 


bq zUUb Greece 4lz-U/-y-Ub-tjr tjU^obb/ 1 


81.4 


91.9 


82.6 


93.1 


82.5 


92.7 


Cr» Alio 1 QQ>1 f^+QT CI l7flQOKn 

rGAUS 1 yo4 oaiy/ bU/UoyOU 


81.3 


93.9 


82.6 


95.1 


82.5 


94.8 


La itaiy iyyb Kvou/yb Arozozu4 


81.4 


93.3 


83.8 


94.5 


83.5 


94.2 


La itaiy zuu i zzy/u i Arozozuo 


80.6 


92 


82.5 


92.9 


82.2 


92.6 


bo inula zUUz Jbo Aroobyi4 


96.7 


96.3 


97.4 


97.5 


97.2 


97.2 


a r> „ i „ _i ; _ onm nQ7n 

Vbo inaia zu iu uy/u 


98.9 


98.8 


100 


100 


99.8 


99.7 


• bo inula z(Jiu bKioo 


98.9 


98.8 


99.8 


99.7 


100 


100 


• bo inula zUUo b/4b 


100 


100 


98.9 


98.8 


98.9 


98.8 


D n inri onnc di id\/*2 cconnc/iQ 
bO mu ZUUb KUdvo brZUUD4y 


99 


98.8 


99.9 


100 


99.7 


99.7 


bq inula zuu/ brvou uuyo 14/ / 


98.1 


97.5 


98.6 


98.8 


98.4 


98.5 


oA 1 1/oU/ 1 A ZUU/ UUOoOD i y 


78.6 


93.3 


79.4 


94.5 


79.2 


94.2 


oi zuu/ oai i/nyb uwooobzu 


78.8 


93.6 


79.6 


94.8 


79.4 


94.5 


ro Argentina iyy4Aioi Loouoo 


77.7 


89.6 


79.4 


90.8 


79.2 


90.5 


ro i naiianu zuud L-ivirUoy ay / u / / oo 


77.8 


90.2 


80.3 


91.4 


80 


91.1 


DO Ur\ ZUUZ L*r- 1 Ar44t5ooZ 


78 


90.8 


80 


92 


79.7 


91.7 


bO ZUU 1 Ur\ rr- 1 Ar4Z/ IZ4 


78 


90.8 


79.2 


92 


79.1 


91.7 


mu Kussia zuuo umsKUo/44z oUo/ / i oy 


77.9 


89.9 


79.7 


90.9 


79.5 


90.6 


Hu Janan 2003 02/92 D86264 


74.6 


86.5 


77.4 


87.7 


77 


of A 


Dog Japan 2004 RAC/DG5 AB526247 


78.5 


89.9 


79.5 


91.1 


79.3 


90.8 


FeAus 1984 Cat2 EU708961 


78.1 


89.9 


79.5 


91.1 


79.3 


90.8 


Hu Germany 2008 GER126/08 FJ747619 


76.2 
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Fig. 3. Phylogenetic analysis of predicted VP7 amino acid sequences of bovine B-46, 0970, and BR-1 33 (G3 strains) with representative 
bovine G3 strains as well as other human and animal G3 strains. The tree was rotoed with Hu India 2009 116E3D (AB081593)RV strain 
sequences. The positions of bovine B-46, 0970, and BR-1 33 are indcated with black circles (•). The percent VP7 nucleotide and amino 
acid identities of the bovine B-46, 0970, and BR-133 strains with representative bovine G3 strains as well as other human and animal 
G3 strains are also presented. 



Phylogenetic analysis 

Phylogenetic analysis of the B-46, 0970, and BR-133 G3 
strains demonstrated that these three bovine strains 
clustered with RUBV3 and J-63 (bovine) and Erv-80 
(equine). This finding indicated the genetic proximity of 
Indian bovine G3 strains to equine RVs (Fig. 3). On the 
other hand, other bovine RVA G3 strains from the UK 
(CP-1, PP-1, and UKg9P) clustered with human, canine, 
feline, and porcine G3 strains (Fig. 3). 

Discussion 

Bovine RVA are associated with high mortality and 



morbidity among young calves, leading to substantial 
economic loss to livestock industry throughout the world, 
particularly in developing countries like India. Previous 
studies from India have shown a predominance of G10 in 
contrast to other G types (G6 or G8) among diarrhoeic 
calves [14,20,22,27,28,34]. However, a report from 
eastern India found a predominance of G6 rather than G10 
as G6 was detected in 89.2% cases and G10 only in 4.6% 
cases [11]. Two new genotypes, G3 (3%) and G15 (2.3%), 
were also discovered in bovine calves in this study. Among 
the P types, P[l 1 ] was predominant (94.6%) followed by 
P[3] (3%), P[2 1 ] ( 1 .5%), and P[ 14] (0.7%). This group also 
detected two unusual combinations, G15P[1 1] (1/130) and 
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G 1 5P[2 1 ] (2/ 1 3 0), in the calves [10]. Prior to these reports, 
Varshney et al. [33] discovered G3 types in diarrhoeic 
calves (10.7%) from central and southern India. The 
frequency of G3 was higher than the G6 or G8 genotypes 
and was second to the predominant G10 type. All these 
studies indicated that G type distribution varies according 
to geographical region. Similar patterns among bovine RV 
G6 populations have been observed in south Ireland [2]. 

G3 RV has been detected in canines, pigs, equines, and 
lapines. Unlike G6, G8, or G10 strains, G3 are rare in 
bovine calves [25]. To date, only two bovine RVA G3 
strains from India (J-63 and RUBV3) and three from the 
UK have been characterized [5,11,25]. In the present study, 
the most prevalent genotypes were G3P[11] and mixed 
genotypes G3+G10P[11], G3G8P[11] and G3G8 P[l] 
P[ll] and G3G10 P[1]P[11], which is uncommon for 
bovines. This indicates a changing pattern of genotype 
distribution from G6, G8, and G10 to a single G3 type or 
G8/G10 in combination with G3 as G3G8 or G3G10 
among Indian bovine populations. Occurrence of the G3 
genotype in central India was also recorded earlier by 
Varshney et al. [33], which further indicates that the G3 
strain is continuously maintained in this region and might 
have spread to other parts of the country given that J-63 and 
RUBV3 are very close in genetic proximity. Variations at 
nucleotide and amino acid levels showed that J-63 and 
B-46 are from the same region but are genetically distinct. 

Observations from the present study are in contrast to 
ones from previous reports [11,14,20,22,27,28,33,34] on 
the predominance of G10 and G6. During the present study 
none of the isolate was typed as G6 while a number of 
isolates were geno typed as G10 type but detected only in 
combination of G3 type. Genotyping failures in RVs has 
been reported earlier due to mismatching/mis-priming 
[1,2,24]; however, all the 33 RV isolates obtained during 
the present study were typable using genotyping PCR 
primers published by Gouvea et al. [12,13], Iturriza- 
Gomara et al. [16] and Taniguchi et al. [31]. The result 
confirms that no other genotype is circulating in the 
regions surveyed during the present study and strains are 
genetically more stable with no variations in the primer 
binding region that normally leads to failures in 
genotyping [1,2,24]. 

The results of present study clearly indicate that G3 in alone 
or combination with G10/G8 is predominant among bovine 
populations in India and is now the most prevalent genotype 
in many parts of the country. No early report from India 
identified combinations of G3 with G8 or Gl 0. This could be 
either due to the use of only G6, G8, and G10 typing primers 
(the predominant G types) or re-assortment/interspecies 
transmission, which may have occurred more recently in 
bovine populations. Atypical detection of G3 or mixed G3 
types in bovine populations indicates a complex interplay 
between the virus, host, and environment. Further studies to 



examine these relationships are warranted. 

Comparison of the predicted amino acid sequences for 
B-46, 0970, and BR- 1 33 G3 RVs to the sequences of other 
VP7 G3 strains showed that Indian bovine G3 strains are 
more conserved in comparison to UK bovine strains. The 
Indian bovine G3 strains appear to represent a 
re-assortment of equine (Erv80) and more recently 
emerged bovine G3 strains. These findings stress the need 
to further investigate any change in genotype distribution 
and the emergence of unusual genotypes through sequence 
analysis. 

In conclusion, detection of G3 or mixed G3G10/G3G8 
types in bovine populations provided significant 
information about the increasing circulation of G3 RV 
genotype in central and northern parts of India. Sequence 
analysis of the G3 strains uncovered evidence of multiple 
species re-assortment events. We observed changing 
patterns of bovine RVA genotypes with increased 
incidences of mixed genotypes containing G3 types. Our 
findings contribute to the growing body of evidence 
suggesting that interspecies transmission and re-assortment 
events between RVs can occur in nature. Therefore, 
examining the distribution of G and P types among bovine 
populations is important for understanding RV ecology 
and the mechanisms by which they evolve, cross the 
species barrier, exchange genes during re-assortment, and 
mutate via the accumulation of single point mutations or 
genetic recombination. Development of novel vaccines to 
effectively combat bovine diarrhoea due to RV will require 
incorporation of the emerging G3 types. 
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